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ABSTRACT 

Palomar 5 (Pal 5) is a faint halo globular cluster associated with narrow tidal tails. It is a useful 
system to understand the process of tidal dissolution, as well as to constrain the potential of the 
Milky Way. A well-determined orbit for Pal 5 would enable detailed study of these open questions. 
We present here the hrst CCD-based proper motion measurement of Pal 5 obtained using SDSS as a 
first epoch and new LET/LBC images as a second, giving a baseline of 15 years. We perform relative 
astrometry, using SDSS as a distortion-free reference, and images of the cluster and also of the Pal 5 
stream for the derivation of the distortion correction for LBC. The reference frame is made up of 
background galaxies. We correct for differential chromatic refraction using relations obtained from 
SDSS colors as well as from flux-calibrated spectra, finding that the correction relations for stars and 
for galaxies are different. We obtain fia = —2.296 ± 0.186 mas/yr and fis = —2.257 ± 0.181 mas/yr 
for the proper motion of Pal 5. We use this motion, and the publicly available code galpy^ to model 
the disruption of Pal 5 in different Milky Way models consisting of a bulge, a disk and a spherical 
dark matter halo. Our fits to the observed stream properties (streak and radial velocity gradient) 
result in a preference for a relatively large Pal 5 distance of around 24 kpc. A slightly larger absolute 
proper motion than what we measure also results in better matches but the best solutions need a 
change in distance. We find that a spherical Milky Way model, with Vq = 220 km/s and V 20 kpc, i-e., 
approximately at the apocenter of Pal 5, of 218 km/s, can match the data well, at least for our choice 
of disk and bulge parametrization. 

Subject headings: Galaxy: fundamental parameters, proper motions, globular clusters: individual: 
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1. INTRODUCTION 

The globul ar clusters of the Milky Way have a ra nge of 
lumin osities (|Harrislll996ll . and Palomar 5 (Pal 5) (lAbelll 
1195511 is one of the faintest. It is located in the halo, at a 
distance of ~ 22 kpc from the sun, an d is associated with 
thin t idal tails of > 22° in length (|Odenkirchen et HI 
120011 iGrillmair fc DionatosI l2006al l. Pal 5 shows un¬ 
ambiguously that tidal disruption is an important pro¬ 
cess for globular clusters, and offers the possibility to 
study tidal disruption in detail at a crucial phase. Other 
globular clusters such as GD-1 appear to have already 
disrupted to thin tails without a clear par ent object 
(jGrillmair fc Dionato^ l2006bl iGrillmai^ 120141 1 . At least 
17% of halo stars show chemic al signs of once havin g been 
members of globular clusters (|Martell et al.l[2Mlll . Thus 
globular clusters are an important contribute r to the stel¬ 
lar halo. According to iDehnen et al.l (1200411 . Pal 5 will 
likely be totally disrupted during its next disk passage. 
However, the orbit of Pal 5 is currently not known well 
enough to strongly constrain its disruption time. 

Streams are useful in determining the shape of 
the dark matter halo because the stars in the tails, 

* Based on LBT data. The LBT is an international collabora¬ 
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LBT Corporation partners are: The Ohio State University, and 
The Research Corporation, on behalf of The University of Notre 
Dame, University of Minnesota and University of Virginia; LBT 
Beteiligungsgesellschaft, Germany, representing the Max-Planck 
Society, the Astrophysical Institute Potsdam, and Heidelberg 
University; the University of Arizona on behalf of the Arizona 
university system; and Istituto Nazionale di Astrofisica, Italy. 

^ Department of Astronomy, University of Virginia, Char¬ 
lottesville, 3530 McCormick Road, VA 22904-4325, USA 
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once stripped from the progenitor, feel negligible in¬ 
fluence from this parent body, and can be used as 
test particles in determining the underlying Milky 
Way p otential (e.g.LTohnston et al.lll999l:lMaiewski et alJ 
200^ iLaw fc Maiewskil I201QI: iKoposov et al.l I201QI: 


^uudt iLaw y iViaiewsKii izuiut ijaoposov et ai.i i^uiut 
Bonaca et al. 1 12014 iPrice- Whelan et al.l l2014f) There¬ 


fore, Pal 5 and its stream offer t he possibility to constrain 
the potential of the Milky Wav (lOd enkirche n et al.ll200.H 

120091: iPearson et al.ll2015l: iKupner et al.ll2015[ l. However, 

existing work on various Milky Way streams offer no 
real conclusion about halo shape. For instance, even the 
well-studied Sagittarius stream has been used to argu e 
for diff erent halo shapes: sp herical (llbata et al.l 1200111 . 
oblate llohnston et al.ll2 nn5li. prolate (iHel mi 2004 1. and 
triaxial ( Law et al .1 1200911 . iLaw fc Maiewskil (|201Clt) fit 
data from the Sagittarius stream, obtaining a slightly 
triaxial Milky Way halo (not outside the expectations 
of A-(-cold dark matter-based structure formation; here¬ 
after ACDM), but with a minor axis that lies within 
the Galactic disk. Such a misalignment is difficult to 
reconcile with galaxy formation models and is expected 
to be ui istable due to torques (iDebattista et al.l I201.H 
see also iVera-Ciro fc Helinl l201ar Modeling of other 
streams usually resu lts in a preference for a spherical or 
slightly oblate halo (iKoposov et al]l2010l : iBowden et al.l 
120151 IlGipper et al.ll2ni,5[ l. Further tests on the shape of 
the halo are thus warranted. 

Streams are also useful tools for mea suring the en¬ 
close d mass of the Milky Way (e.g. iGibbons et al.l 
[201j). The M ilky Way mass is still poorly constrained. 
IGibbons et al.l (120141 1 obtain a value of 6 x 10^^ Mq us¬ 
ing the Sgr stream, nearly a factor of three less than the 


































































2 


FRITZ & KALLIVAYALIL 


mass obtained bv iBovla n-Kolch in et alJ (I2013D : 16 x 10^^ 
Mq using Leo I; see lWang et al.l (|2015D for a recent com¬ 
pilation of Milky Way mass measurements. Due to the 
fact that it is relatively close-by, Pal 5 can also be used 
to constrain the Milky Way disk. 

FinallVi Pal 5 show s tantalizing gaps in its stream 
(|Carlberg et al.l l2012tl. These may be caused by sub- 
structure (|Dehnen et al.l 1200411 . like giant molecula r 
clouds, spiral arms or dark subhalos (ICarlbergI 1201211 . 
or by the epicyc lic nature of t idal d i sruption even in 
a smooth ha l o (iKupper et al.l 120081 lJust et al. l 


Kiipper et al 


Kiipper et al. 


j20^ IMastrobuono-Battisti et all 


2009; 


2012 


1201511 . Gaps are well detectable only in 


thin (cold) st reams, which makes Pal 5 useful to st udy 
substructure (|Yoon et al.ll201lUCarlberg et al.ll2012ll . 

To shed light on all of these topics a proper motion 
measurement would be useful. A proper motion in com¬ 
bination with the known radial velocity of the cluster 
gives Pal 5’s present-day kinetic energy directly, which 
enables better orbit constraints, and in turn, fewer free 
parameters in modeling the tidal disruption of the sys¬ 
tem. Currently, there are thre e photographic plate -based 
proper motion measurements (jDinescu et al.l[l99^ which 
contradict each other. In this work, we obtain the first 
CCD-based proper motion of the Pal 5 cluster, using two 
large field-of-view imagers, SDSS and the Large Binoc¬ 
ular Camera. In Section [2] we present the data. We 
describe the measurement of the proper motion of Pal 5 
in Section |3l In Section |4] we use the obtained proper 
motion to explore the Milky Way potential and the or¬ 
bit and distance of Pal 5. We discuss and conclude in 
Section [S] 

2. DATA SET 

In this section we describe the data used to measure 
the proper motion of Pal 5. 


2.1. SDSS 

We use SDSS data release DR9 (|Ahn et al.l 1201211 as 
our first epoch of data. Subsequent data releases do not 
contain new imaging. We retrieve all clean objects in 
the Pal 5 stream area from the Photo table. This in¬ 
cludes generous margins around the stream debris lo¬ 
cation, and also covers our second epoch imaging area. 
Selecting clean objects means that we exclude dupli¬ 
cates, sources with obvious deblending and/or interpo¬ 
lation problems, spurious detections, and stars close to 
the edge. The clean criterion does not work perfectly; 
in addition we later exclude sources using other crite¬ 
ria, as discussed in Section 13751 In the case of stars we 
use PSF magnitudes and in the case of galaxies we use 
cmodel magnitudes. SDSS uses cmodel magnitudes to 
separate stars from galaxies. Stars (ProbPSFband =1) 
are defined as PSFmagband“cmodelmagband< 0.145. For 
the DCR correction we use the measured magnitudes, 
see Section o For the purpose of selecting of Pal 5 
stars fSection 13.31) . we use extinction-corrected magni- 
tudes, relying on the SDSS-provided correction based on 
iSchlegel et al.l (jl998l l. We use the given SDSS positions 
and their errors as a base. We then iterate on these po¬ 
sitions to better fit the differential chromatic refraction 
(DCR) corr ecti on a s we ll as the error model, as discussed 
in Sections 13.21 and 13.41 The epoch of the SDSS data is 
1999, with an average MJD= 51269. 


2.2. LBC/LBT 

For the second epoch o f imaging we use the Large 
Binocular Camera (LBC) (jCiallongo et al.l 120081 1 at the 
Large Binocular Telescope (LBT). These data were ob¬ 
tained on 2014 July 1, giving a 15-year baseline between 
the first and second epoch of imaging. All LBC images 
have exposure times of 40 seconds and cover a field of 
view of about 23' x 25'. The field is divided into four 
detectors, three upright at the bottom and one prone at 
the top. Each individual detector covers 7.8' x 17.6'. We 
obtained 12 images of the cluster and 41 images of the 
stream. The images of the cluster overlap generously, 
while the images of the stream (usually two images per 
position) have only small overlaps between them, see Fig- 
ure[Il We use the r-band data of the red eye of the LBT. 
The r-band is the reddest of the high SNR LBC bands 
and is thus less influenced by atmospheric effects such as 
DCR. The pixel scale of LBC is 0.226". The typically 
achieved FWHM is 0.9". 

We use standard data reduction techniques on the LBC 
data: we obtained skyflats, biases, and from the combi¬ 
nation of the two, bad-pixel masks. We do not correct 
for cosmic rays or other artifacts (such as those around 
very bright stars) in the individual images. Instead we 
use outlier rejection to exclude bad pixels in our source 
lists. Since we have between 2 and 12 images of each 
pointing, this is robust and does not reduce the number 
of sources in a relevant way. Stars close to very bright 
stars are usually already excluded from the clean SDSS 
source list. Thus, for our current purposes there is little 
point in applying more advanced data reduction tech¬ 
niques to be able to extract more sources close to very 
bright stars from the LBC data. We do not combine 
the individual LBC images into stacks, as position er¬ 
rors can be better estimated by using each of the indi¬ 
vidual images. Furthermore, image-stacking can often 
degrade the astrometry, when the images cover different 
pointings (especially when the pointing difference is only 
slightly smaller than the field of view of the instrument, 
see e.g. iCillessen et al.l[2(10^ . and whe n the pixel scale 
varie s significantly over the field of view (|Giallongo et al.l 
[2noMlBellini fc Bedinllmiill . 

3. DERIVING THE PROPER MOTION OF THE PAL 5 
CLUSTER 

It is in principle possible to measure the proper motion 
of the stream from our data. However, due to the small 
contrast of the stream compared to the background (see 
e.g. Figure [T] top) more aerial coverage is necessary than 
has been obtained here, and this will instead be the fo¬ 
cus of subsequent work. In this paper we focus on the 
measurement of the cluster proper motion. 

3.1. Obtaining pixel positions 

Since we use SDSS as our first epoch, we are con¬ 
strained to use only SDSS-detected sources in our anal¬ 
ysis. Starting with these as an input list, we use a pre¬ 
liminary version of our coordinate transformations (Sec¬ 
tion [T4]) to obtain the expected pixel positions of SDSS 
sources in our LBC images. We then search within a 
radius of 4 pixels (~ 0.9”) for the local maximunjj. 

^ We use here for many basic steps: dpuser, see 

http: //www.mpe.mpg.de/~ott/dpuser/index.html 




















































Palomar 5 


3 



R.A. [°] 



R.A. [“] 


Fig. 1. — Map of LBC coverage of Pal 5. Top, the stream area: 
the blue dots show the stars which fall in the Pal 5 matched filter, 
and thus are likely members of the Pal 5 system. The area of the 
dots is proportional to the weight within the filter. The cluster 
Pal 5 is around 229.02/-0.12. The red polygons show the field of 
view of one LBC pointing. We cover the area with 22 different 
pointings. Bottom, the cluster itself: blue dots show stars in the 
matched filter, the central hole is caused by source blending in 
SDSS. The green circle (r= 0.13°) marks the outer boundary of 
the cluster as used in this work. The red rectangles show the four 
detectors of the four different pointings on the cluster. The orange 
dots show the reference galaxies. Their area is antiproportional to 
their position uncertainty. 


This search box is sufficient to find all objects except 
for high proper motion stars, which are foreground stars 
and not the targets of the present analysis. We then 
fit a two-dimensional Gaussian to the 4-pixel half-width 
box around that maximum, leaving both widths and the 
orientation of the major axis of the Gaussian as free pa¬ 
rameters. We exclude very small and large widths, which 
likely result from failed fits. At this point we do not dis¬ 
criminate between stars and galaxies. We use two vari¬ 
ants of the LBC images, the standard reduced images, 
and images which we obtain by smoothing these with a 
Gaussian of FWHM= 4 pixels. There are no systematic 
differences in the positions obtained from these two kinds 
of images. However, we use the positions obtained from 
the smoothed images, since in this case the fits fail more 
rarely. The fits fail only for ^ 6 % of the sources, which 
are usually faint. The smoothing also reduces the dif¬ 
ference in shape between our simple source models and 
the actual sources that in reality are more complex. It 


is unlikely that our source shape model is limiting our 
accuracy. Our error floo r fSection 13.31) is consistent with 
the error floor of SDSS (|Pier et alJl2003H . Also the t otal 
error is dominated by the SDSS errors, see Section I.S.bl 
which shows that our simple shape model is sufhcient to 
obtain positions. 


3.2. Differential ehromatic refraetion correction 

The atmosphere of the earth deflects light with an an¬ 
gle of 

a = a'tan((j) ( 1 ) 

in which (j is the angle from zenith and a' follows from 
the index of refraction, n(A): 


, _ njXr - 1 

2n(A)2 


( 2 ) 


a is 57” for ( = 45°. Most of the deflection, a, is cor¬ 
rected for automatic allv in any linear transformation, 
see Section 13.41 and iFritz et all (|2010fl . However, the 
strength of the refraction also depends on the wavelength 
of the lig ht. While that de pendence is small in the near- 
infrar ed (iFritz et al.ll201C )fl. the effect is larger in the op¬ 
tical (IKaczmarczik et al.ll2009H . In the optical: 


[n(A)-l]-10® = 64.328-f 


29498.1 
146- (1/A)2 


255.4 

41- (1/A)2 


(3) 


is valid (jFiliDpenkol Il982t) . While the absolute refrac¬ 
tion does not matter here, DGR causes a wavelength- 
dependent angular offset, /3offset, between objects of dif¬ 
ferent wavelengths. To correct for /3offset, it is neces¬ 
sary to know the effective wavelength (Agff) of an object 
within a given ba ndpass. For photon counting detectors 
the d efinition is (|Fritz et al.l 120111 : iTokunaga fc Vacc^ 

[IMl): 


fA^F^(A)S(A) 

fAF^{A)S{A) 


(4) 


Therein S{A) is the full transmission curve of the 
filter used, including the atmospheric transmissior 0 , 
and F\{A) is a flux-calibrated spectrum. Thus, flux- 
calibrated spectra are necessary for the calculation of 
effective wavelengths. From these spectra we derive for 
each object class the color-Agff relation. These relations 
can be converted into color-/3offset relations for a given 
(, se e Figure [2] SPSS uses flux- calibrated s tar sp ectra 
from [Gunn &: Strvk^ (|1983ll . see IPier et al.l (12001^ . We 
also use these spectra for our stars to derive the color-Agff 
relation. To construct a color for our sources, we use the 
f-band. The Tband is the best band to derive a color-Agff 
relation since it is closest in log(A) to the r-band. It also 
contains less strong features than the 5 -band. 

We next test whether galaxies follow the same color- 
Agff relation as the stars. Since SDSS spectra do not cover 
the full imaging magnitude range (down to r = 21.5), we 
instead use spect ra from the zCOSMOS bright sampl^l 
(jLillv et al.l I2009H to obtain Agff (and both datasets in 
the calculation of colors, as described below). Compared 


® see http://classic.sdss.org/dr7/instruments/imager/ for SDSS. 
® Based on observations made with ESO Telescopes at the La 
Silla Parana! Observatory under program ID 175.A-0839. 
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to some other public samples this sample has the advan¬ 
tage that the spectra also cover the Aband. The two 
data sets (SDSS and zCOSMOS) differ in their treat¬ 
ment of Galactic extinction. The Galactic extinction is 
calibrated out in the stellar spectra while it is present in 
the galaxy spectra. In principle, it would be best to use 
extinction-free spectra in both cases, and to use these to 
also calculate the effect of extinction. However, we deem 
that this is unnecessary in our case: the extinction to 
Pal 5 {Ar—Ai = 0.026 to 0.048) and towards the zCOS¬ 
MOS field {Ar—Ai = 0.011 to 0.016) is small compared 
to the color range of stars {r — i ^ —0.4 to 1) and galaxies 
{r — i ^ 0.1 to 1.2). Furthermore, the extinction-related 
DCR correction follows approximately the same function 
of color as the intrinsic color-/3offset correction. 

For stars we derive a tight color-/3offset relation, see 
Figure [2] (we exclude some noisy very late-ty pe dwarfs). 
Howe ver, a linear relation, as used by SDSS (|Pier et al.l 
1200311 . does not fit the full color range. We instead use 
a spline relation. For normal-colored stars the spline is 
close to the linear relation used by SDSS. We set /^offset 
of stars with r — i = 0 to zero. 

For the galaxy relation we use 303 galaxies from the 
zCOSMOS sample to calculate Aeff. This is a sufficient 
number to not be limited by statistics. The high S/N 
part of these spectra lies redward of 0.562 fim. To ob¬ 
tain an estimate of the spectra blueward of this out to 
0.538 fj,m (the blue edge of the r-band filter), we fit the 
range between 0.562 and 0.577 fj,m. and extrapolate. We 
note that this is only a small extrapolation, as this range 
contains only a small fraction of the total r-band flux. 

We then quantitatively examine the resulting (r — i)- 
Aeff relation for the galaxies, as a function of both the 
SDSS colors and the synthetic colors calculated from the 
spectra. The (r — i)-Aeff relation has more scatter and 
contains more outliers when we use the SDSS colors. 
That is to be expected, since the SDSS data is nois¬ 
ier. We therefore do not use SDSS colors of individual 
galaxies for the determination of the (r — i)-Aeff rela¬ 
tion. The median color difference between the synthetic 
and the SDSS galaxy colors is (r — i)syn — {'^ — *)sdss = 
0.092 ± 0.009. This offset could be physical, because 
SDSS measures the full galaxy color while the zCOS¬ 
MOS slits cover only the (usually redder) galaxy centers. 
However, the offset is more likely a result of a spectral 
calibration issue. We explore the astrometric effects of 
both possibilities below. 

We exclude 17 galaxies from our sample, due to noise or 
extreme emission lines. We examine the resulting (r — i)- 
PoSset relations in Figured] The (r —i)-/3offset relation for 
galaxies is noisier than for stars. This is expected from 
the fact that galaxies have more variable extinction. The 
galaxy (r — i) — /Joffset relation is also offset from that of 
stars. This offset is caused mainly by the fact that galax¬ 
ies consist of many types of stars. To test this we tested 
how exti nction-free, s i ngle s tellar populations, such as 
those of iCoelho et al.l (1200711 . look on such a diagram, 
and found that a single-population sequence is itself off¬ 
set from our stellar sequence by about the same amount 
as the galaxy spectra. 

We fit the galaxy {r — i) — /^offset relation with a linear 
model. A fourth order polynomial reduces the scatter 
by only 6 %, and we therefore stay with a simple linear 
model. We set the few very red {r — i > 1.1) and blue 



r-l (AB) 


Fig. 2. — Differential chromatic refraction as function oi r — i 
color. /3offset is the shift to zenith for = 45°. The shift for stars 
with r — 2 = 0 is set to zero. The black line is the relation used by 
SDSS. We obtain for stars a similar relation (violet line). However, 
our relation for galaxies (orange line) is clearly offset from this 
relation. 

{r — i < 0.1) sources to have a constant offset given by 
the value at the border of their defining range, and fit the 
remaining galaxy data with a linear model, obtaining 

PoSset = (11.17±0.24)-(35.76±1.00)x(r-i) [mas] (5) 

for ^ = 45°. The difference between model and data has 
a scatter of 3.9 mas. 

We then test how a spectral calibration issue would 
affect the astrometry. In this case we assume that the 
correct synthetic flux, Fcai, can be obtained in the fol¬ 
lowing way from input spectra, Flnput: 

Feed = -Finput X A"^ (6) 

From fitting the above equation with the given flux in¬ 
puts, we get X = —0.5. For the color term in Equation 5 
we use the median for galaxies in the area of our Pal 5 
data, which is r — i = 0.48. For this median color, /^offset 
changes by only 0.3 mas when we use Fcai instead of 
Finput from Equation [6] This is comparable to the l-cr 
error of Equation O The change is small because the 
changes in r — i and /3offset work to cancel each other 
out. Because the difference is not significant, it does not 
matter for the astrometric accuracy whether the spectra 
are offset due to a calibration issue or by a physical ef¬ 
fect. The fact that our DGR zero-point is anchored by 
stars and not an ‘average galaxy’ has the consequence 
that the positions of our galaxies are offset on average. 
However, we are not interested in perfect absolute posi¬ 
tions for galaxies, we only want a reference system which 
does not move due to DGR, which we achieve with the 
correction from Equation [5] 

For the application of the correction we first divide 
the sources into stars and galaxies as in Section 13.31 
Note that quasars are stars according to this classifica¬ 
tion. Since t heir spectra are differen t from both stars 
and galaxies (|Kaczmarczik et ^l2009f ) our DCR correc¬ 
tion will likely not work on them. Thus, we do not use 
them for tests and ignore them in this work. We also 
need to use a different DCR treatment for SDSS versus 
LHC positions. The SDSS positions are already in equa¬ 
torial coordinates and are corrected for DCR. Since we 
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have derived here a better DCR correction we subtract 
the applied DCR correction from the SDSS positions and 
then apply our new one. In the case of LBC, we derive 
a linear coordinate transformation from pixel to equato¬ 
rial coordinates for each image and detector, to convert 
(doffset from mas to pixels. Our transformation assumes 
perfect alignment of the pixel axes with R.A. and Dec. 
While this is not exactly true, the deviations from these 
assumptions are smaller than 2% in total. We then ap¬ 
ply the pixel DCR corrections to the pixel positions. The 
corrected pixel positions are then converted back to equa¬ 
torial coordinates (see Section IXTl) . 

The mean difference between the raw and the DCR- 
corrected motions is Sfia = —0.03 mas/yr and 6^s = 0.38 
mas/yr. From the 3.9 mas scatter, and factoring in the 
number of galaxies and the time baseline, we derive an er¬ 
ror of less than 0.02 mas/yr for the final DCR correction. 
Even if this error is underestimated due to unknown sys- 
tematics it is not significant compared to our total proper 
motion error of about 0.18 mas/yr per dimension. 

3.3. Reference frame and target star selection 

To obtain proper motions it is necessary to identify 
member stars and to measure their motions relative to 
a reference system with a known velocity. In principle, 
point sources are preferable as a reference since the po¬ 
sition of sourc es with small FW HM can be measured 
more precisely (iFritz et al.ll201Clf) . However, there are no 
good point sources in our case. The density of quasars 
is too low for sufficient precision. The average veloc¬ 
ity of many faint halo stars is larger than our measure¬ 
ment accuracy (see e.g. iBond et al.l (1201011 '). Thus, stars 
can only be used when their velocities are known with 
high precision. There is no reliab le option in our case: 
the Hipparcos astrometry catalog (iHog et al.ll2000ll only 
contains stars that are saturated in our data and their 
density is also too low. P hotographic plate-b ased survey 
catalogs such as those of iMunn et al.l (|2004ll are deeper. 
However, while they have good per-star precision, it is 
possible they contain systematic erro rs. It is risky to use 
a galaxy model (see the example of iPrvor et al] (120101 1 
for the Sagittarius dwarf) to predict the mean velocity 
of the foreground star population, since Pal 5 is distant 
from objects with well measured proper motions like Sgr 
A*. We thus use galaxies as our reference frame: they are 
abundant and we also have better control in estimating 
the measurement errors. 

As a starting point for the galaxy sample we use the 
well vetted morphological classification of SDSS in the 
r-band. A potential worry for the reference galaxy sam¬ 
ple is pollution from stars. We reduce the likelihood of 
this in several ways. The likelihood of pollution is espe¬ 
cially high towards the dense center of Pal 5. Therefore 
we check galaxies within r — 0.07° of the cluster center 
by eye in the SDSS and LBC data. We find that galax¬ 
ies which fall approximately along the cluster CMD se¬ 
quence look star-like and have a star-like FWHM in the 
LBC images. We therefore exclude all galaxies in this 
CMD region, as well as any other galaxies which appear 
star-like in shape or FWHM, from our reference sample. 

Furthermore, we use the LBC images to refine our 
selection. We use SEx t ractor togeth er with PSFex 
((Bertin fc Arnouti 119961 iBertinI 1201111 in the three- 
step process outlined by iKoposov et al.l (|2015ll . We 


use as the criteri o n for galaxies I SPREAD_M0DEL I > 
0.003 (iDesai et HI 120121 : lAnnunziatella et all 120131 
IKoposov et al.ll20T5l) . Only 4% of the galaxies selected 
according to the previous criteria are rejected by this 
criterion. We verify by eye that most of these appear 
star-like in the LBC images. Since all of these criteria 
need to be fulfilled at once, our galaxy selection is more 
robust than when using, say, only the SDSS morpholog¬ 
ical classification. 

Next we use variants of these criteria to estimate the 
uncertainty caused by the galaxy selection. The SDSS 
criterion is PSFmag-cmodelmag= 5mag = 0.145 for each 
band separately. We test this criterion in three differ¬ 
ent ways: first, we use the information from all bands 
with enough signal {g and f-band) in addition to the 
r-band, which is our main workhorse; second, we add a 
significance criterion compared to the magnitude error of 
galaxies (tr cmodelmag): dgai = {dma-g) / {o' cmodelmag); 
we use (Tgai = 5, which means we only use objects which 
are significantly brighter when fit with a galaxy model 
than when fit with a stellar PSE model. The LBC selec¬ 
tion remains the same as described above. In both cases 
we calculate the change to the final proper motion by 
the selection, we obtain changes of order 5 /Xgai ~ |0.04| 
mas/yr per dimension. That is small compared to the 
final total accuracy of 0.18 mas/yr per dimension (Sec¬ 
tion [331). 

To identify stars belon ging to Pal 5 , we use the 
matched filter technique (jRockosi et al.l 1200211 . The 
matched filter is defined in the following way: 

W = E(C'MD)targetobject/S(C'MD)background (7) 

Yj{CMD) x above is the density of stars in bins of color 
and magnitude in a CMD, either the CMD of the target 
(Pal 5) or the background. For the CMD we use the r 
and g-bands, which have the highest SNR. We first select 
all stars in SDSS within 8.4’ o f the cluster center (the 
tidal radius of Pal 5 according to lRockosi et al.ll2002ll , see 
Figure [3] We then exclude stars that do not match the 
majority of Pal 5 stars in the color magnitude diagram 
via outlier rejection. 

There are several options in constructing the Pal 5 
matched filter, W. One option is to do simple binning, as 
in a Hess diagram. However, that is not optimal, because 
it is influenced by shot noise in the case of a small num¬ 
ber of stars. Another option, binning with smoothing, 
makes the shot noise smaller; however, it also makes the 
cluster sequences wider than they are in reality. Instead 
we model the cluster CMD as much as possible: we fit the 
main sequence, the (sub)-giant branch and the horizontal 
branch with different Gaussians whose widths, heights, 
and positions vary smoothly as a function of magnitude. 
For the blue straggler region Gaussian modeling does not 
work well. There we use a smoothed Hess diagram in¬ 
stead. 

For the background we choose a 4 square degree region 
that is offset from the Pal 5 cluster, with borders of 225° 
and 227° in R.A. and 1.9° and 3.9° in Dec. This region is 
large enough that shot noise is not important and small 
enough that the Milky Way background is still close to 
identical to the Pal 5 region. We smooth the Hess di¬ 
agram of the background slightly to further reduce the 
influence of shot noise. The final filter W is shown in 
Figure [3] There, and in general, we use a magnitude 
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Fig. 3.— CMD of Pal 5 region, using extinction-corrected SDSS 
PSFmag values, and excluding very red sources {g — t > 1.2). 
The blue dots show Pal 5 stars. The gray background shows the 
matched filter W which we use to identify Pal 5 stars. The red 
contour levels in linear steps show the density in the background 
field. 


cut-off of r = 22 which is approximately the complete¬ 
ness limit of SDSS for point sources. The precise cut-off 
is not important since stars close to this limit have large 
motion errors such that they are given negligible weight 
in the subsequent analysis. 

To get membership probabilities, p, for each star, first 
a relative weight (w) is obtained by checking in the 
matrix W the value at the magnitude and color of a 
given star. We then construct the background density 
(E(ry)back) and cluster density (E(r(;)pai 5 ) from these 
weights. For the cluster, we measure the density ra¬ 
dially in bins from the center, see Figure IH E(ui)pai 5 
dips in the cluster center because the high source den¬ 
sity reduces the number of clean SDSS sources. We 
mask out this central depression and fit the density of 
the other radial bins with a polynomial. To measure 
the density of the background we use a 2° box, centered 
on the cluster, from which we mask out the cluster and 
the stream. The radial membership fraction is defined 
as /pai5(?’) = (E(r(;)pai5 - E(r(;)back)/E(w;)pai5. We then 
obtain a probability p for every star using the follow¬ 
ing procedure: for each /pai 5 (r') we solve the following 
equation for niiim: 


/pal5(?’) 



( 8 ) 


We then draw for each star a random number (between 
0 and 1) and multiply it by its w. When that product is 
larger than wum the star is used. This way we create 10 
different random samples, in order to test whether the 
uncertainties in membership are important. We do not 
use stars outside of r = 13'. At 13' / = 45%. 

The error contribution to the proper motion from this 
Pal 5 star selection, as calculated from the scatter of the 
10 different samples, is small. We obtain a contribution 
of only 0.046 mas/yr. This is less than the contribution 
from distortion and image registration (see Section 13.41) . 
We did not test whether galaxies misclassified as stars 
are a relevant problem. However, because we select Pal 5 


Ppal 

0%<Ppal 5<10% 





0 5 10 15 20 


r ['] 

Fig. 4.— Top: CMD of Pal 5 (r<13’). The probability for mem¬ 
bership depends on the position in the CMD (the value of the 
matched filter W there, see Figure and the distance from the 
cluster center. Stars with zero probability are not plotted. Bottom: 
Radial density profile of Pal 5 after application of the weighted fil¬ 
ter. The central depression is caused by source confusion. We fit 
the other data with a polynomial. 


stars in the CMD and because there are more Pal 5 stars 
than galaxies in the cluster region, this effect should be 
small. 

3.4. Distortion correction and image registration 

Distortion is oft en the main error source in astrometry 
(|Fritz et al.ll2010ri . There are two main ways to estimate 
the distortion. Firstly, one can use many dithered obser¬ 
vations to infer the distortion field of a given instrument 
without external referen c e fsee e.g. p Vnderso n fc Kin^ 
l2003t iTripoe et al.l 120081 iBellini fc BedinI I2010D . Sec- 
ondly, one can use a distortion-free source list of the ob¬ 
served area. While the first strategy has higher position 
precision (a position error of crAi n = 10 mas in the case 
of LBC, see lBellini fc BedinI 1201011 it is necessary to be 
able to apply it to all epochs. This is a disadvantage for 
us because we use different instruments. While such a 
correction would be possible for the LBC, it is more dif¬ 
ficult in the case of SDSS. It is further unclear whether 
improvements can be m ade over the exis ting SDSS dis¬ 
tortion calibration from iPier et al.l (|2003f) . 

For our correctio n we use the SDSS astrometry as refer¬ 
ence. According to iPier et ahl (120031) the SDSS internal 
precision is about 25 mas. Thus, we expect a similar 
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accuracy. In addition to the distortion, it is also neces¬ 
sary to solve for the image registration (linear terms). In 
our definition the image registration contains the linear 
parameters and the equatorial coordinate of the central 
pixel, which changes according to the pointing. Unlike 
the distortion terms which can be stabl e over time, the 
linea r parameters are usually not stable (|Bellini fc BedinI 
I20IOH . due to different effects, for example, the changing 
airmass. Constant i mage scale parameters limit the ac¬ 
curacy to about 1" (jBellini fc BedinI r20ir)H . Thus, it is 
necessary to solve for all image registration parameters 
in every image. 

In practice we use the following process: first, we ap¬ 
ply the distortion correction to the raw pixels to obtain 
undistorted pixel coordinate^ Second, we apply the 
image registration to the corrected pixels to obtain sky 
coordinates. For the distortion correction we use the fol¬ 
lowing equation, wherein x' and y' are measured relative 
to pixel (1049, 2304) respectively, the approximate center 
of the detectors: 

Xcon = X + ai x'^ + 02 x'y' + 03 y'^ + 04 x'^ 

-1-05 x '^^ y ' + 06 x ' y ''^ + 07 y '^ 

2/cor = y + bix''^+ 62 x'y' + h 2/'^ + &4 x'^ 

+&5 x'^y' + be x'2/'^ -I- 67 y'^ 

Using higher order corrections than this do not improve 
the residuals. Our residual scatter for bright galaxies is 
close to the internal precision of SDSS. The distortion 
correction does not contain linear terms because these 
are contained in our image registration transformation: 

R.A. — Cl -j- C2 Xcor F C3 year 

Dec. = di+d 2 Xcor + ds 2/cor (10) 

We test three different variants of pixel coordinate 
transformation: firstly, time-dependent (separate for 
each image) and detector-dependent distortion param¬ 
eters and image registration which, as stated, vary 
with time and detector. Secondly, time-independent 
but detector-dependent distortion and image-registration 
parameters. Thirdly, time-independent but detector- 
dependent distortion parameters and time-dependent im¬ 
age registration, which are coupled between the different 
detectors. Specifically, in the third case, we first fit the 
same time-dependent image registration parameters to 
all detectors. We then fit a constant skew and shift term 
for each detector except for detector one which we use as 
the “master” detector. This case has clearly larger resid- 
uals than the other tw o cases, confirming the finding of 
iBellini fc BedinI (|201Clll that using all LBC detectors to¬ 
gether leads to larger errors than when the detectors are 
treated separately. In a variant of case three we couple 
the offsets between the detectors to the variable image 
scale. This assumes that only the airmass changes the 
angular distances between the detectors. In this variant 
the residuals are reduced but are still larger than when 

^ Before that we exchange x and y definitions for detector 4. 
Then they are approximately aligned with the others detectors. 


the detectors are treated independently. Thus the detec¬ 
tor distances depend on at least one other variable, pos¬ 
sibly temperature variations. We do not use case three 
because of its larger errors. 

The first two cases give consistent results; that means 
the distortion is stable over the 80 minute time-scale 
of our observations. We decide to use case two: inde¬ 
pendent detectors with a constant distortion solution for 
each, but time-variable image registration terms. Com¬ 
pared to case one it has the advantage that it has fewer 
free parameters. 

The small number of galaxies in the very core of Pal 5 
(see Figured]) affects the precision of the linear parame¬ 
ters adversely. We therefore use the cluster stars as addi¬ 
tional reference sources. Since the cluster stars move, we 
add the displacement in R.A. and Dec as additional pa¬ 
rameters. These two parameters can be easily converted 
to the proper motion of Pal 5, by dividing through the 
time baseline between the two epochs. This is easily pos¬ 
sible because the time baseline is effectively the same for 
all stars, since there is only a one day difference in the 
observing epochs of the Pal 5 cluster in SDSS, and only 
a few minute difference in LBC. 

Adding the distortion terms and the image registration 
terms together we fit 334 parameters for each of the four 
detectors. The number of sources used in the fit depends 
on the detector and the drawing of Pal 5 stars. There 
are about 11,970 sources per detector. There are always 
more than enough so urces to fit the pa rameters. For the 
fitting we use mpf it ()Markwardtll20n9l l . For the position 
errors of the objects we use the SDSS errors (ctXidsdss, 
called raErr and decErr in SDSS) as a base. From that 
we calculate the full error: 

crXiD = \/(ctAidsdss X factY + add? (11) 

add is the systematic, magnitude-independent error, 
which is not included in the ctXsdss id of DR9. fact 
can stand for LBT-related errors and/or serve to correct 
o'Asdssid- As a first guess we use add = 20 mas and 
fact =1. In our first derivation of the transformation 
terms we fit only galaxies. We calculate the residual {R) 
of the positions relative to the fit, and exclude outliers 

with R = ^R^ + Rf^ > 10 <7. A possible explanation 

for outliers are objects misclassified in SDSS, bad pixels, 
and in the case of complex galaxy morphologies, center¬ 
ing problems. 

We then analyze the residuals (after outlier rejection) 
in the following way: we bin these galaxies according to 
(tXidsdss in 20 equally-populated bins, see Figure[SJ In 
each bin we determine the scatter of the residuals be¬ 
tween the two epochs. For robust outlier-rejection we do 
not use the standard deviation, but the median devia¬ 
tion. We scale it by a factor of 1.483, to make its value 
identical with the standard deviation under the assump¬ 
tion of a Gaussian distribution. We fit the error model 
(Equation [TT|) to the points, finding that it is a good 
representation of our data as can be seen in Figure [SJ 
Only the point for the largest errors is clearly offset, pos¬ 
sibly due to the fact that the error for the point may be 
overestimated. We exclude it; due to its large error, its 
influence is very small. 

We repeat the fit several times with iterative outlier 
rejection. In the process we also add the Pal 5 stars 
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Fig. 5.— Residual position scatter of LBC sources compared to 
SDSS sources. We use 1.483X the median deviation as a robust 
scatter measure. The line shows our error model. 


to the data sample. We use for them the same error 
model, since the error should not depend on the nature 
of the source. The parameters of our final error model 
are add =17.21 and fact = 1.186, see Figure[5j It seems 
slightly conservative especially in add for the objects in 
the Pal 5 cluster center field. However, it is difficult to 
measure add reliably from this relatively small area. For 
the final outlier rejection we choose R = 3.21. With this 
cut the final x^/d.o./. = 1, which implies that our sample 
does not contain further outliers. The distribution is not 
fully Gaussian; there are a few more stars with very small 
errors and relatively large residuals than expected from a 
Gaussian distribution. As a test we also model the errors 
by fitting only cluster stars or only galaxies, finding that 
the obtained motions do depend slightly on the error 
model used. From the scatter in these different results, 
we obtain a motion error of 0.04 mas/yr, which we fold 
into our final errors below. 

Overall we have enough coverage (Figure [IJ to ob¬ 
tain a robust distortion solution, see Figure [51 The dis- 
tortion is qualitatively simi lar to the one obtai ned by 
iBellini fc BedinI (l2010ll and iFabrizio et al.l (l2014f l. The 
differences can be explained by the use of a different fil¬ 
ter and camera and by time variability on the timescale 
of years. Our solution has a ID accuracy of 17 mas over 
the ful l area. (Using only t he cluster the accuracy is 13 
mas.) IFabrizio et al.l (I2014D obtain ctAid = 6.7 mas as 
precision for LBG. The difference is likely caused by the 
fact that our error estimate contains also th e SDSS er¬ 
rors, which alone are the size of our uXid ([Pier et al.l 

Hooi. 


3.5. Proper motion of Pal 5 

To obtain the final motion we combine the Pal 5 dis¬ 
placements of the four different detectors weighted by 
their errors. For the final errors, we first correct for the 
fact that the displacements obtained from the different 
LBC images are correlated. The reason is that they all 
use the same SDSS positions, which dominate the errors. 
We first multiply the errors by where N = 10.5, is 
the average number of detected galaxies. This assumes 
that only SDSS contributes to the error budget. To in¬ 
clude non-zero LBC errors, we need to multiply this by 


another factor, M < 1. From the strength of the posi¬ 
tion correlations between the different LBC images we 
derive M = 0.944. Combining the two factors, we derive 
ap,a = 0.176 mas and apig = 0.172 mas from the galaxy 
and star errors (uXid). To this error term we also add 
the uncertainties from image distortion and registration 
('Section lHlTl) . and the unc erta inty due to the Pal 5 mem¬ 
bership selection fSection lT3)) . Other error sources con¬ 
tribute negligibly. Considering all these effects we obtain 
fia = —2.296 ± 0.186 mas/yr and pLg = —2.257 ± 0.181 
mas/yr. 

In Figure 13 we compare our measurement (the first to 
be obtained from CCD data) with the v alues from photo¬ 
graphic plate studies in the literature (iSchweitzer et al.l 
119931 : iScholz et al.lll998L K. M . Cudworth 1998, unpub- 
lished. see lDinescu et al 111999( 1. The scatter between the 
photographic plate-based motions is larger than their 
quoted errors. The mean of these measurements is 
roughly consistent with our value. 

4. CONSTRAINTS ON THE MILKY WAY HALO 

We now use the 6D phase space properties of Pal 5, 
i.e., its position, distance, radial velocity, and measured 
proper mot ion, and the publicly available code galpy 
(lBovvll2?)T5h . to model the disruption of the Pal 5 sys¬ 
tem. We compare the predictions of the model to the 
measured location of the stream debris and the observed 
radial velocity gradient along the stream. 


4.1. Used parameters 

While the position of Pal 5 is kno wn to sufficien t preci¬ 
sion, its distance is more uncertain. iHarrisI (|1996ll derives 
a distan ce of 23.2 kpc froiu the h o rizontal branch, as is 
used by lOdenkirchen et aki (l2009tl. iVivas fc ZinnI (I2006f) 
derive 22.3 kpc from RRLyrae. iDotter et al.l (|2011|) ob- 
tain 20.9 kpc0 from fitting an HST CMD. It is unclear 
from the respective error analyses whether these values 
are in agreement with each other. In order to explore a 
realistic range, we usually use the smallest (20.9 kpc) and 
largest (23.2 kpc) values as a bound. As our “standard” 
value we use the mean of these values: 22.05 kpc. 

The radial velocity o f the cluster, W, is -58.7 
km/s (lOdenkirchen et aki 120021) . A comparison with 
iKuzma et al.l (|2015| ) shows that the error is approxi¬ 
mately 0.8 km/s. The measured dispersion of the cluster 
is (J^aib(rad) = 1.1 ±0.2 km/s, which upon correction for 
inflati on caused by binary stars is CTpai 5 (rad) ~ 0.3±0.15 
km/s (lOdenkirchen et al1l2002D . Due to mass loss, the 
dispe rsion of Pal 5 was higher in the past ([Dehnen et al.l 
120(141) . Therefore, we use a value of (Tpai 5 (rad) = 0.4 
km/s in our modeling below. 

For the stream we use two observables, the stream 
streak on the sky and the radial velocity gradient along 
the stream. The radia l velocity gradient along the 
stream was measure d bv lOdenkirchen et al. (120091) and 
IKuzma et al.l (|2015D . lOdenkirchen et al.l ( 20091) mea¬ 
sured velocities for Pal 5 stars using high resolution spec¬ 
troscopy. In their work they use only spectroscopically 
confirmed giants. Four giants are clearly not associated 
with Pal 5 because of their radial velocities. In the case 


® Note that the distance moduli cited in IDotter et al.l 1120111) are 
not corrected for extinction (private communication), and therefore 
need to be corrected before use. 
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TABLE 1 



STREAM 

STREAK 

POSITIONS 


R.A. [°] 

Dec. [° 

o-Dec. [°] 


241.48 

6.41 

0.09 


240.98 

6.15 

0.09 


240.48 

6.20 

0.09 


239.98 

5.81 

0.09 


239.48 

5.64 

0.09 


238.48 

5.38 

0.09 


237.98 

5.14 

0.09 


233.61 

3.17 

0.06 


233.11 

2.88 

0.06 


232.61 

2.54 

0.06 

1998 

232.11 

2.23 

0.06 


231.61 

2.04 

0.06 


231.11 

1.56 

0.06 


230.11 

0.85 

0.06 


229.61 

0.54 

0.06 


228.48 

-0.77 

0.11 


228.11 

-1.16 

0.14 


227.73 

-1.28 

0.11 


227.23 

-2.03 

0.17 


226.55 

-2.59 

0.14 


Fig. 7.— Proper motion measurements of Pal 5. 


of two others it is less clear whether they are members, 
binaries whose velocity is influenced by their partner, or 
a re not members. 

iKuzma et alJ ()2015f l use (mostly) medium resolution 
spectroscopy to obtain giants over a larger range of the 
stream. For selecting the Pal 5 stars they use spectro¬ 
scopic giant-dwarf and metallicity discriminators. Fi¬ 
nally, they exclude stars above -30 km/s and below -70 
km/s. The stars at the largest distances from the cluster 
are so mewhat off the stream locatio n (as measured by 
us and iGrillmair fc DionatosI (j2006ari . While this maybe 
worrisome, their simulation shows that pollution by non¬ 
stream stars should not be a relevant problem. There ¬ 
fore , we use all the memb ers of lOdenkirchen et al.l (l2009f l 
and iKuzma et al.l (j2015H together with the cluster veloc¬ 
ity to obtain the gradient as a function of R.A. Since 
we are interested in the mean velocity of the stream, 
we add in quadrature to the individual me asurement er¬ 
rors the stream dispersion of 2.1 km/s (jKuzma et al.l 
[Ml). First, we also try to incl u de th e two dubious 
members from lOdenkirchen et ^ (I2009D . but that re¬ 
sults in y^/d.o.f. = 116.2/56. Excluding them, we ob¬ 
tain x^/d.o.f. = 47.4/54. Thus, the velocities of these 
two stars are very likely influenced by some non-stream 
motion (such as binarity). We therefore exclude them 
and obtain a gradient of 1.21 ± 0.09 kms“^deg~^ in R.A. 
The gradient is identical within the e rrors when using all 
stars, or o nly the stars f rom e ither lOdenkirchen et ^ 
([Ml or IKuzma et al.l (1201511 . Most stars at large 
distances from the cluster are on the same side of 
the stream, which causes some additional uncertainty. 
To account for this as well as other potential uncer¬ 
tainties, we enlarge the errors slightly and use 1.20 ± 
0.15 kms“^deg~^ as gradient when we compare models 
with the gradient. 

As for the stream streak, because it is difficult to esti¬ 
mate the errors from published maps, we construct our 


own. We retrieve all stars from SDSS for the stream 
and adjacent areas. We begin by binning the stars in 
quadratic bins of 3.75’ in length. For map one we give 
all the stars the same weight. For map two we weight 
the stars with the filter constructed in Section 13.31 To 
account for varying surface density along the stream, we 
divide map two by map one. We treat bins with less than 
20 stars as bad pixels and interpolate over them. To in¬ 
crease the signal to noise ratio of this map we smooth 
it with a 30’ Gaussian, obtaining the final map of Fig¬ 
ure [H We sample the map at 30’ R.A. intervals, and fit 
a Gaussian plus a gradient at each location. As a first- 
guess error for each position, we ensure that the relative 
errors account for the variable density of the stream. We 
treat the leading and trailing streams separately, since 
the stream follows a complex shape (an ‘S’) close to the 
cluster. A linear model is sufficient for the short leading 
stream, while we use a quadratic model for the longer 
trailing stream. Finally, we rescale the previous errors by 
a common factor to achieve x^/d.o./. = 1 for the leading 
and trailing streams separately. We present the obtained 
positions in Table [TJ Our pos itions agree well with the 
recently published positions of lKlipoer et al.l (j2015H 
The velocity and position of the sun are also necessary 
inputs. For the distance Ro of the sun to the Galac¬ 
tic center w e combine thr e e rece nt high-accuracy mea¬ 
surements: iDekanv et al.l (1201311 obtained 8. 33 ± 0.15 
from RR Lyrae in the bulge. iReid et al.l (|2014ll obtained 
8.34 ± 0.14 kpc from parallaxes and velocities of masers 
in a su bstantial part of the disk, and lGhatzopoulos et al.l 
(l2015fl obtained 8.27 ± 0.13 kpc from a nuclear cluster 
model fit to radi al velocities and proper motions from 
iFritz et ^ (|2014ll . The stated errors combine the sta¬ 
tistical and systematic errors. Combining these consis¬ 
tent measurements weighted by their errors we obtain 
Rq = 8.31 ± 0.08 kpc. Th is Rn is also consis tent with 
most older measurements ([Genzel et al.ll2010h . but has 
a smaller error. We also use the following value for the 
distance z of the sun r elative to the mid plane, which is 
z = 0.02 ± 0.007 kpc (Oo^ 120071 : iMaiaess et al.ll200l : 
iBuckner fc Froebrichll2014ll . 
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Fig. 8.— Matched filter map of the Pal 5 stream. The red points 
show the stream positions which we use in our modeling. 


TABLE 2 

Parameters of Pal 5 cluster, stream and Sun 


parameter 


value 


R.A. 

Dec. 

distance 

MaPalS 

M(5Pal5 

yatream (rad) gradient (R.A.) 

rsun 

ZSun 

Vsun(r) 

VSun(0) 

Vsun( 2 ) 


229.018° 

-0.124° 

20.9 to 23.2 kpc 
—2.296 it 0.186 mas/yr 
-2.257 ±0.181 mas/yr 
-58.7 ±0.8 km/s 
0.4 km/s 

1.2 ± 0.15kms“^deg“^ 
8.31 ± 0.08 kpc 
0.02 ± 0.007 kpc 
-11.1 ±0.7 km/s 

251.3 ±2.6 km/s 

7.3 ± 0.4 km/s 


It is standard to dissect the solar motion into the cir¬ 
cular rotation around the Galactic Center. Vn, and its 
peculiar motion, UsR l|Kerr &: Lvnden-Bel]|ll986D . Since 
Vo an d Vr.su are both uncertain at possibly the 20 km/s 
leveljSchdnrich et al.ll201(il : iBovv et 311120121 : iReid et id] 

I2014D . a better approach for our purposes is to use the 
motion of the sun more directly. The proper motion of 
the Sun is visib le as reflex motion in th e proper mo¬ 
tion of Sgr A* (IReid &: Brunthaleil 120041) . Combining 
the proper motion with the distance gives Vsun(</') and 
Vsun(-2)- This Vsun (z) is Consistent, but less precise, 
than the Vsun (z) determined from the LSR. We there¬ 
fore u se the Rsun(0) from above, and the lSchdnrich et ahl 
(IMIl values for Vsun (z) and Vsun (r ). These solar 
motion components are not disputed ( Schonrich et al.l 
l201Ct IBovv h all 120121) . Our final solar values are: 
Vsun(c (I), z) = -11.1 ± 0.7,251.3 ± 2.6, 7.3 ± 0.4. (Due 
to our coordinate system definition V(r) is reversed in 
comparison to U in the LSR.) Because of our derivation 
method, Vsun( 0 ) and i?o have a strong positive corre¬ 
lation. These and all other parameters for the Pal 5 
system (apart from the stream streak) are summarized 
in Table [2j 


4.2. Galaxy model setup 


Thanks to the mass-ratios involved, msat << Afhost, 
satellite debris streams provide strong constraints on 
the host potential, and for computational simplicity, 
streams are often modeled as generally following or¬ 
bits. Iloweyer, streams do not strictly fo l low orbits, see 
e.e. iJ ohnstoiil (119981): iHelmi et al.1 (119991): IDehnen et al.1 

e ; KiiDDer et al.l (I2012D : iSanders fc BinnevI (I2013D : 

(I2014D . N-body simulations have shown that 
stream-orbit misalignments can be significant, and in¬ 
crease with the mass of the cluster a nd the eccen¬ 
tricity of the orbit (jLiuc et al.l 120131 ). While N- 
body simulations are a conceptually sim ple method to 
account for stream-orbit mis alignment (|Dehnen et al.l 
12004 iLaw fc Maiewskil I2010D . they are computation¬ 
ally expensive, making it difficult to explore parame¬ 
ter space. There are now several alternative theoret¬ 
ical approache s to obta i ning st ream position s and ve- 

who use 


locities, e.g. ISandersI (120141 ) 
an action-angle f o rmalism, and 
iLane et ahl (120121): iKiipper et al 


Bov 


L position ! 


Varghese et al.l (I 2011 D: 

(I 2012 D: iGibbons et al l 


( 20141) : iPearson et al.1 ( ^151) : iKiiDper et al.l (I2015D . who 

use a combination of a restricted three-body approach 
and N-body simulations. Here , we use the publicly avail¬ 
able action-angle formalism of iBovvl (j2014D , galpy (jBovvl 

nnn). 

For testing purposes, we use the following Milky Way 
models. Our main model (Standard) i s very simi¬ 
lar to the model used by IPearson et al.l (j2015D . and 
to a common parametrizati on of the Galactic poten¬ 
tial (jAllen fc Santillaiil 119911), which uses the lAU rec - 
ommendation of 220 km/s (jKerr fc Lvnden-Belllll98^ . 
It has three components: bulge, disk and halo. For the 
bulge we use a (spherical) iHernaui^ (119901) parametriza - 
tion; for the disk we use a iMivamoto fc Nagail (jl975f ) 
parametrization; for the halo we use a spherical loga¬ 
rithmic parametrization. This and the parameters of the 
three components are close to identical with the spher¬ 
ical halo model of IPearson et al.l (I2015D . We list them 
in Table |3l Besides this main model we use two vari¬ 
ants of it: Small and Massive. In these variants the 
mass of bulge and disk are identical to what is in the 
Standard model, but we vary the ma s s of t he halo, sim¬ 
ilar to the approach of lKupper et al.l (|2015D . We choose 
the mass of the halo such that the rotation curve in¬ 
side of 30 kpc is in rough agreement with other obser¬ 
vations at the solar radius (lMcM illan| [ 20 Tl|: IBovv et al.l 
l 20 I 2 l:IPenarrubia et al.llMl4lR.eid et al.ll20I 4D and out to 


about 30 kpc (iSakamoto et al.l 20031: Koposov et al.ll2?)Tfll : 


iKafle et 1111120121 : IKiipper et al. 1^151) . We~ ignore con¬ 

straints further out, because Pal 5 does not enter these 
distances. 

A logarithmic potential is probably not a good 
parametrizati on of the full halo of the Milky Way. A n 
NFW-profile (|Navarro et ahllTo^ IKiipper et al.ll2015D is 
probably a better representation of the Milky Way halo. 
However, the difference between a logarithmic and an 
NFW halo is small over the limited radial range probed 
by Pal 5. The chosen models have Vq in the range of 
210 km/s to 235 km/s. A better measure to characterize 
the models in o ur context is the velo city at the apocen- 
ter of Pal 5, see IKiipper et al.l (I2015D . For simplicity we 
assume here that the apocenter is at 20 kpc in the disk 
plane. The velocity there (V 20 ) is between 192 km/s and 
249 km/s, while our Standard model has a value of 218 
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r [kpc] 


Fig. 9. — The rotation curves of the three potentials. All use the 
same disk and bulge, but the (spherical) halo mass is allowed to 
vary. 

km/s. The rotation curves of the three models are shown 
in Figure m 

4.3. Comparison with models 

After initializing with the current-day 6D measure¬ 
ments of the Pal 5 system, we model the disruption of the 
cluster using galpy, an orbit-integration-based approxi¬ 
mation to calcul ate ac t ion-an gle coordinates described 
in Appendix A of iBovvl (|2014ll . We produce both a lead¬ 
ing and a trailing stream, assuming 5 Gyr since the start 
of disruption (the default in galpy). Given the relatively 
short length of the observed stream 22°), which is 
probably mainly limited by SDSS coverage, we do not 
attempt to constrain the disruption time. However, in 
our tests we find that the disruption time has no influ¬ 
ence on the specific stream observables used below. 

We also do not attempt to fit the radial velocity dis¬ 
persion of the progenitor, apa\ 5 {rad). crpai 5 {rad) has the 
largest influence on the offset of the stream relative to 
the cluster at points closest to the cluster center. The 
influence of crpais {rad) on the stream at further distances 
from the cluster center is relatively small. Therefore, 
we keep the disruption time and the dispersion fixed 
and explore other areas of parameter space: the halo 
mass, which is connected with the model parametriza- 
tion: Standard, Massive or Small, Rq, the velocity of 
the Sun in the Galactic plane (Vsun(</<)), the distance 
of Pal 5, and the proper motion. (The position and 
radial velocity of the cluster and the solar parameters 
tsun, Vsun(?’), and Vsun(t) are well-known in compari¬ 
son, and their errors affect our an alysis here negligibly. ) 
One finding in our w ork and that of iPearson et al.l (1201511 , 
IKupper et al.l (120151 1 is that when the proper motions in 
the two dimensions are changed by different multiples of 
their errors, the resulting stream has a different shape 
than what is observed. Therefore, when we change the 
proper motion, we increase (or decrease) the proper mo¬ 
tion in both dimensions by a common multiple of their 
errors. 

We find, as expected, that all these parameters make a 
difference, but three of them are clearly more important 
than the rest: the halo mass, the distance of Pal 5, and 
the proper motion. 



E.A, [•] 


Fig. 10. — Comparison of the observed stream streak with differ¬ 
ent models. The crosses mark the location (o, <5) of the observed 
stream debris (Figure [8ll. The lines show the different models. 
When not indicated otherwise they use the standard parameters, 
Vq = 220 km/s, 22.05 kpc, and the standard proper motion (pm) 
of -2.30/-2.26 mas/yr. In the two cases shown, both dimensions 
of the proper motion have been changed equally in the same di¬ 
rection. The details of the model parameters are explained in the 
text. 

In Figure [TU] we compare the stream streak with the 
output of the stream disruption models. Specifically, the 
interpolated track location is what is plotted here as the 
output of the stream disruption model. Due to its shorter 
length, the leading stream is much less constraining than 
the trailing stream. As discussed, the explored solar pa¬ 
rameters are not important for the stream streak: when 
we decrease Rq and Vsun(?i>) both by 1 cr, the streak 
moves up only by 0.4° at the very end of the stream 
(R.A.= 241°). The Standard model with Vq =220 km/s 
and V 20 =192 km/s produces a trailing stream which 
clearly has lower Declination values than the observed 
one. A smaller V 20 J as in the Small model, decreases the 
discrepancy with the data. However, our Small model 
is still clearly offset from the stream. Even smaller V 20 
values are unlikely based on constrains from other data 
(iSakamoto et aLM^llO.II: iKoposov et al.l 120101: iKafie et ^ 
l2012HKuDDer et al.ll2015l) . While we cannot exclude that 
other changes in the details of the potential setup, such 
as flattening, can result in a good fit, we now explore the 
impact of the phase space parameters of Pal 5 itself: its 
proper motion and its distance. 

We find that larger (more negative) proper motions 
and/or larger values of distance obtain better consistency 
with the stream streak, see Figure [TUI Both imply that 
a larger velocity for Pal 5 in physical units fits the data 
better. Good agreement is achieved for either a distance 
of 24.2 kpc or a change in the proper motion of about 
— 1.5(7. Because the distance to the cluster has no for¬ 
mal uncertainty, it is difficult to assess the significance of 
changes in the distance. Gurrently, both the lower bound 
and the upper bound of the distance we use are equally 
likely. Therefore, given that the mean value of our mea¬ 
sured proper motion requires larger distances in order to 
fit the stream better, we interpret this as evidence for a 
larger Pal 5 distance. The fact that the Small models 
fit the data better can be interpreted as a preference for 
a small V 20 and thus a small halo mass. However, since 
other parameters, i.e., distance and proper motion, have 
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TABLE 3 

Model parameters 


Model name 

Component 

fiiVo) 

Vio 

parameters 



Standard 

Bulge 

0.307 

121.9 

c = 0.7 kpc 



Standard 

Disk 

0.496 

154.9 

a = 6.5 kpc, 

b = 

0.26 kpc 

Standard 

Halo 

0.197 

97.65 

ry = 12 kpc 



Massive 

Bulge 

0.270 

121.9 

c = 0.7 kpc 



Massive 

Disk 

0.4345 

154.9 

a = 6.5 kpc, 

b = 

0.26 kpc 

Massive 

Halo 

0.2955 

127.75 

ry = 12 kpc 



Small 

Bulge 

0.337 

121.9 

c = 0.7 kpc 



Small 

Disk 

0.555 

154.9 

a = 6.5 kpc. 

b = 

0.26 kpc 

Small 

Halo 

0.119 

72.44 

ry = 12 kpc 





R.A. 

Fig. 11.— We s h ow al l velocities of the stream members from 
lUdenkirchen et al.l (120091) (violet) and IIKuzma et al.l [201511 (ma- 
genta). For the latter, we do not show the targets that overlap with 
IQdenkirchen et al.l (|2Q09I ). The errors include a contribution from 
the measurement errors as well as the dispersion of the stream. We 
calculate the gradient as function of R.A. from these data points; 
black line with hatched 1 a error area. The colored lines show the 
different models. When not indicated they use the standard pa¬ 
rameters of Vb = 220 km/s, 22.05 kpc and -2.30/-2.26 mas/yr. 
For the two labeled cases, both dimensions of the proper motion 
have been changed equally in the same direction. The details are 
explained in the text. 


more leverage within their uncertainties, we do not want 
to overstate this preference, especially in light of adding 
further stream constraints to the modeling, as discussed 
below. 

The observed radial velocity gradient along the stream 
provides an additional constraint. In Figure [TT] we show 
the corresponding agreement in this space. Again, solar 
parameters are not important, a reduction by 1 cr reduces 
the velocity only by 3 km/s at R.A.= 241°. The impor¬ 
tant parameters are found to be similar to the stream 
streak case, i.e., cases which fit the stream streak well 
also fit the radial velocity gradient well. 

Again either a larger distance, or larger (more nega¬ 
tive) proper motion are preferred. The best fitting val¬ 
ues this time are similar as in the stream streak case, 
about -1.3(7 for the proper motion and 24.2 kpc for the 
distance. Thus, the Standard model can fit all the data 
with a distance of 24.2 kpc. This is on the high end 
of the currently measured distance, and since the pub¬ 
lished distance values do not report errors, it is difficult 
to make a more quantitative statement here. In the case 
of changing the proper motion the fit is somewhat worse. 


because a different proper motion is necessary to fit the 
radial velocities versus the stream streak, while a single 
distance can fit both equally well. 

In Figure [10] and [Tl] it is also visible that the Small 
model can fit both data sets better for the standard dis¬ 
tance, although it is a worse fit than changing either 
proper motion or distance values. To quantify this we 
calculate the of the fit to both the stream streak and 
the radial velocity gradient for several parameter choices. 
The is calculated as follows. For the stream streak 
we calculate, for each data point, the smallest distance 
to each model streak and use the distance to calculate 
the for each data point using our error in the stream 
position. Similarly, for the radial velocity, we calculate 
for each star the smallest distance to the model gradi¬ 
ent to get the for each star. The radial velocity error 
consists of the measurement error of each star and the 
stream dispersion of 2.1 km/s added in quadrature. 

In Figure 112] we show the total x^ as a function of dis¬ 
tance for our three Milky Way models. The smallest 
value (of 70/74) is that of the Standard model at 24.2 
kpc. The x^ values for the other models are higher be¬ 
cause in these cases different distances are necessary to 
minimize the for the stream streak versus that of the 
radial velocity. This effect causes the Small model to 
fit worse than the Standard model. When varying the 
distance, the Small model fares better than the Massive 
model. The reverse is true when we vary the proper mo¬ 
tion, even to very large values. Thus, both Small and 
Massive models are approximately equivalent, and lead 
to worse fits than the Standard model. As one addi¬ 
tional test on the dependence of the details of the Galaxy 
m odel we use th e inter polated centerline best fit values 
of IKiioper et al.l (|2015ll (Column 4 of Table 6 ) 0 . This 
work uses a similar disk and bulge parametrization as 
used here, but an NFW halo. They marginalize over 
proper motion and distance, obtaining best-fit values of 
23.12 kpc and -2.50/-2.44 mas/yr. When using exactly 
these parameters (NFW halo, distance and proper mo¬ 
tion), we obtain the lowest y^ of all investigated models, 
of 65.8/72, even when we add to this y^ another term 
that takes into account the difference between the model 
proper motion and our measured proper motion, which 
results in a total y^ of 68/74. To remind the reader, this 
is slightly better than the y^ of Standard model with 
24.2 kpc. 

® We choose that model over the other models of this work 
because in galpy only spherical NFW halos are currently imple¬ 
mented. 
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. Massive 235 
Standard 220 
Small 210 


distance [kpc] 


Fig. 12. — shown as a function of distance, for our three 
models. All models use our measured proper motion of -2.30/-2.26 
mas/yr. The calculation of the uses both the distance from the 
stream streak and from the stream radial velocities. 


Overall we find that none of the tested models can 
fit the data, even approximately well, for the standard 
distance of 22.05 kpc and proper motion. Varying both 
quantities, and within the limitations of our Milky Way 
model search, we find two Milky Way models that fit the 
data well: our Standard model which has Vn = 220 km/s 
and ¥ 20 = 192 km/s, and one model of iKfipper et al.l 
(|2015f) . which has Vq = 238 km/s and V 20 = 236 km/s. 
A more thorough Milky Way model search would be eas¬ 
ier if a better distance to Pal 5, with an error estimate, 
were available. 

The constraints on some orbit parameters are better: 
the three best fitting cases, large distance, large abso¬ 
lute proper motion, bot h with Standard mo del, and in¬ 
terpolated centerline of iKhpoer et al.l (|2ni5h with their 
best-fit parameters, all indicate approximately consistent 
pericenters between 7.74 and 8.28 kpc, respectively. In 
contrast, we obtain looser constraints on the apocen- 
ter because it mainly depends on the current distance 
of Pal 5, which is itself close to apocenter. We obtain 
between 17.76 and 19.86 kpc for the apocenters of these 
three cases, respectively. 

5. DISCUSSION AND CONCLUSIONS 

We have measured the proper motion of Pal 5, using 
SDSS DR9 data and LBT/LBC data, separated by ^ 15 
years, obtaining a proper motion of jia = —2.296 ±0.186 
mas/yr and fia = —2.257 ± 0.181 mas/yr. We discuss 
this first in the context of the disrupti on of Pal 5. We ob¬ 
tain a pericenter distance of ~ 8 kpc. lOdenkirchen et al.l 
(j2003li determin e the pericent e r of t he next passage to 
be about 6 kpc. iDehnen et al.l (12004( 1 use a pericenter of 

5.5 kpc in their N-body simulations to predict the fate 
of Pal 5 during the next orbit, obtaining a passage very 
close to the disk plane. Independent of the precise mass 
and shape of the cluster, they obtain that Pal 5 will be 
disru pted in the next disk passage. In their recent mod¬ 
eling (KuppeZir^ ()2015[) obtain a pericenter of about 

7.5 kpc. Therefore, the pericenter of the s tream orbit is 
proba bly larger than in the simulation of IDehnen et al.l 
(|2004ll . Thus, tidal effects and shocks will probably be 
weaker in the next passage, which makes it possible that 
the cluster will not be disrupted in the next passage. 


Now we compare with other wor ks on the Galactic po¬ 
tential. lOdenkirchen et ahl (|2003ll predicted the proper 
motion of Pal 5 using the radial velocity of the clus¬ 
ter and their stream coverage, which includes the cen¬ 
tral 10° of the stream. They use a distance of 23.2 
kpc and a logarithmic (spherical) halo potential (with¬ 
out a disk or bulge), and with a constant Vq of 220 
km/s. From orbit fitting to the stream they obtain 
^.a/ = —2.03/ — 2.045 mas / yr, jus t outside of our 1-cr 
errors. iGrillmair fc DionatosI ((2006a[ ) predict the proper 
motion by using the radial velocity of the cluster, the 
stream streak, Rq = 8.5 kpc, a distance of 23.2 kpc an d 
the Vo =220 km/s potential of I Allen fc SantillanI (il991f) . 
which consists of a bulge, disk and spherical halo, very 
similar to our Standard potential. They obtain a mo¬ 
tion of /iaZ/ra = —2.27/ — 2.19 mas/yr, within our 1 cr 
interval. Both works assume that the stream follows the 
orbit. 

We also comp a re o ur results with those of 
lOdenkirchen et ahl (l2009f) , who compare their ana¬ 
lytic results to a prediction of t he stream-orbi t offse t 
from a full N-body simulation by IDehnen et al.1 (|2004l) . 
They find that this offset becomes important especially 
when the radial velocity gradient is also used in the 
fitting. In that respect, Pal 5 is different from GD-1, 
for which stream-fitting (ignoring stream-orbit offsets) 
seems to do as well as a full N-body model, However, 
that is possibly caused by the potent ia l parametriza- 
tion used by iKonosov et al.l ()201(lli : iBowden et al.1 
because whether a stream follows an orbit is 
a very strong fun c tion o f the global potential model. 
lOdenkirchen et ahl (l2009f l obtain a good match to both, 
the stream streak and the radial velocities, for Vq ~ 220 
km/s. They prefer small values of Vq, citing 200 km/s or 
180 km/s as a better fit to their methods. Based on our 
tests here, it is possible that the fact that their potential 
is spherical explains w hy they prefer a rather small Vq. 

iPearson et al.l (l2015f) use a combination of a restricted 
3-body approach (the streak-line method) to explore pa¬ 
rameter space, followed by more detailed N-body simu¬ 
lations, to study the disruption of Pal 5. They leave the 
proper motion of Pal 5 as a free parameter, and seek to 
reproduce both the stream contours and the radial ve¬ 
locity gradient, given either a spherical Milky Way halo 
model (similar to ours i n Sec tion 14.21) . or the triaxial 
halo of iLaw fc Maiew^ ((20101) . In the triaxial case they 
find a best-fit proper motion of = —5.0/ — 3.7 

mas/yr, inconsistent with our measurement. The best 
fit proper motions for the spherical halo, however, of 
Ma/ltiS = —2.35/— 2.35 mas/yr, are very cons istent with 
ours. The conclusion of IPearson et al.1 (I2015D . is that a 
spherical halo gives a better fit to t he stream contours 
and t he radial velocity gradi ent of lOdenkirchen et al.1 
((20091) than the triaxial halo of iLaw fc Maiewskil ((201(1lb 
with the triaxial halo causing the stream to “fan out” 
much more than observed. Our measured proper mo¬ 
tion is consistent with their predicted motion in such a 
spherical halo for the Milky Way, as is our own model¬ 
ing of the disruption of the cluster. However, our results 
do not rule out the presence of a radial gradient in the 
halo, wherein it changes from more spherical in the in- 
ner parts to more triaxial in the outer parts ( see also 
iLaw fc Maiewskil2010HVera-Ciro fc Helniil[2013i) . 

Very recently. iKiioper et al.l ( 20151 ) used most Pal 5 
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observables (except for proper motion), including the 
stream streak, and the r adial velocity gradient data of 
lOdenkirchen et ^ (l2009f) to fit the parameters of the 
Milky Way halo, the Sun and Pal 5. They used the 
streak-line method, obtaining solar parameters which are 
very similar to the solar parameters used by us. The 
model setup is also similar to ours except that the halo 
is pa rametrized with an NFW potential (|Navarro et al.l 
I1997I1 . and the flattening of the halo is left as a free 
parameter. They obtain a proper motion of ^a/pi-s = 
—2.39to;i7/ — 2.36lQ;i[5 ruas/yr consistent with our ob¬ 
servations. They also obtain a high distance estimate 
of 23.58;°;?^ kpc. As a normalization, they obtain 
K) = 233t}2-7 km/s. The halo is consistent with being 

spherical with = 0.95lQ;i[2- 

Since our proper motion is slig htly smaller than the 
prediction of iKiioper et al.l (|2015h . it will be interesting 
to use this as a prior to obtain new constraints on Vq 
and the dist ance. In genera l, our c onclusions are similar 
to those of iKiioper et al.l (|2ni5h . like the preference 
for a relatively larger distance of Pal 5. However, this 
similarity may be caused by the fact that the bulge 
and disk in our works are nearly identical, which has a 
significant effect in the radial range of Pal 5 (^7—20 
kpc). In particular, both of our studies employ a rather 
massive and extended disk in c o mparison to recen t 
works, see e.g. iKlvpin et al.l (120021) : Rix fc Bovvl (j2013l l: 
iMcMillaiil ()2011ll and I Bow fc Rixl ( 2013D . Such a disk 
not only influences the fla ttening of the potential (see 
e.g. iKoposov et al.l 1201011 but also the enclosed mass 
within the range probed by Pal 5. A further issue of 
both studies is that there is currently no good distance 
measurement for Palomar 5, thus adding another nearly 


free parameter. A good measurement of the distance 
of Pal 5 with an error would be valuable. We will 
explore the effects of a smaller disk and different param- 
eterizations of the halo, the statistical uncertainties in 
the solar parameters, and the proper motion and dis¬ 
tance of Pal 5, together with a full disruption model, to 
investigate the disruption of Pal 5 in a subsequent paper. 
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